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ABSTRACT 



Context. A short observing run with the spectrometer HARPS at the ESO 3.6-m tele- 
' scope was conducted in order to continue exploring the asteroseismic properties of F 

type stars. In fact, Doppler observations of F type on the main sequence are demanding 
and remain currently limited to a single case (HD 49933). Comparison with photomet- 
ric results obtained with the CoRoT mission on similar stars will be possible with an 
enhanced set of observations. 

Aims. We selected the 4th magnitude F8V star HD 203608, in order to investigate the 



oscillating properties of a low-metallicity star of the old galactic disk. 
Methods. A 5-night asteroseismic observation program has been conducted in August 
2006 with HARPS. Spectra were reduced with the on-line data reduction software pro- 
vided by the instrument. A new statistical approach has been developed for extracting 
the significant peaks in the Fourier domain. 

Results. The oscillation spectrum shows a significant excess power in the frequency 
range [1.5, 3.0 mHz]. It exhibits a large spacing about 120.4 /iHz at 2.5 mHz. Variations 
of the large spacing with frequency are clearly identified, which require an adapted 
asymptotic development. The modes identification is based on the unambiguous sig- 
nature of 15 modes with £ = and 1. 

Conclusions. This observation shows the potential diagnostic of asteroseismic con- 
straints. Including them in the stellar modeling enhances significantly the precision 
on the physical parameters of HD 203608, resulting in a much more precise position in 
the HR diagram. The age of the star is now determined in the range 7.25 ± 0.07 Gyr. 
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1. Introduction 



New stable spectrometers dedicated to very precise radial velocity measurements have per- 
mitted rapid progress in observing solar-like oscillations in solar like stars (see for example 
Bedding & Kjclds en 2006| for a recent review). The questions raised with the observation 
of the CoRoT target HD49933 (jMosser et al. 2005|) . an active F5V star, led us to continue 
the observations of F stars. We therefore planned to measure, identify and characterize the 
solar-like oscillations of another low-metallicity F type dwarf star, with a 5-night run. 

The selected target HD 203608 (HIP 105858, HR 8181, 7 Pav) is a F8V star which 
belongs to the group of the Vega-like stars with an infrared excess attributed to the 
presence of circumstellar dust warmed by the central star. Its age is estimated to range 
between 6.5 to 14.5 Gyr ( |Bryden et al. 2006[ ), according to previous work including 
Lachaume et al. (1999)| based on different estimators: comparison to theoretical isochrones, 



rotational velocity, strength of chromospheric calcium emission lines, stellar metallicity, and 
space velocity. Its metallicity [Fe/H]=— 0.65 dex designates it as an old dwarf of the optical 



thick disk according to Edvardsson et al. (1993) but Bensby & Feltzing (2006) describes 
it as a member of the thin disk. 



Pijpers (2003) has already identified HD 203608 as a possible asteroseismic target. Its 



right ascension and declination are favorable for a single-site campaign at La Silla during 
winter. Moreover this star has a low projected rotational velocity (v sini = 2.4±0.5kms~ 1 , 
IReiners fc Schmitt 2003|) which is excellent for Doppler asteroseismology, especially for a 
F star. On the other hand, with broader lines than G stars, F stars are demanding targets 
for ground-based seismic Doppler observations. As a result, HD 203608 is an exciting star 
to compare with a similar star, the active F5V star HD 49933, a main target of the satellite 
C0R0T already studied by spectrovelocimetry (jMosser et al. 2 005). 

Section [2] reports the current status on the physical parameters of HD 203608. 
Observations are presented in Section [31 with the analysis of the time series and of the 
activity signal. The seismic analysis, based on the unambiguous detection of £ = and 
1 modes is exposed in Section 2J Asymptotic parameters are extracted from the Fourier 
spectrum, then individual eigenfrequencies and amplitudes. The modeling of HD 203608 is 
presented in Sect.0 Section [5] is devoted to conclusions. 



2. Stellar parameters 

2.1. Temperature and luminosity 

The atmospheric parameters of HD 203608 have been discussed in several studies 
(see |Cayrel de Strobel et al.T997| . It results a range of 5929 K < T off < 6139 K for 
the effective temperature. Using photometric data, Thevenin (1998) and more recently 
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Jonsell et al. (2005) and del Peloso et al. (2005) adopted a similar value: T e g — 6070 K 



with an uncertainty of ±100 K. But this star is a Vega-like one which exhibits infrared 
excess ( |Mannings fe Barlow 1998] ). The flux absorption and reemission in the infrared is 
difficult to estimate, and its effective temperature may be underestimated. 

The absolute bolometric magnitude was determined using the apparent V magnitude 
(TableQJ combined with its Hipparcos parallax (it — 107.98±0.19 mas. Ivan Leeuwen 2007P . 
The bolometric correction used is from Flower (1996) The resulting luminosity is L/L Q — 
1.39 ±0.13 (Table©. 

2.2. Abundances 

The atmospheric stellar abundance published in the recent literature converges to the value 
—0.67 dex for the iron element with a detailed analysis under LTE assumption. Using NLTE 



computations, Thevenin & Idiart (1999) (hereafter IT99) have proposed to correct both 
this iron LTE abundance and the spectroscopic surface gravity. They have shown that the 
more metal-poor the star, the more important NLTE effects: the ionization balance of iron 
adopted to determine the spectroscopic log g is in error. IT99 having adopted T c s — 6072 
K. Like |Jonsell et al. (2005)] we can use the NLTE values proposed by IT99 for [Fe/H]: 



—0.57 dex. The logc/ value is deduced from the Hipparcos distance. This parameter has a 
large uncertainty (~0.15dex) due mainly to the bolometric correction, and therefore cannot 



be used to constrain efficiently the stellar fundamental parameters. Korn et al. (2003) found 
also NLTE effects on the lines of the iron element in metal-poor stellar atmospheres but in 
a lesser quantity since they adopted for the NLTE computations a higher efficiency of the 
inelastic collisions with hydrogen. In the absence of a realistic theory of inelastic collisions 
with H, we shall adopt a NLTE abundance of —0.60 dex for Fe. In most of metal-poor 
stars the alpha elements present an enhancement compared to iron element abundance, 



up to 0.4 dex. Idiart & Thevenin (2000) have computed NLTE abundances for Ca and Mg 



elements in HD 203608. For that star, Ca appears to follow Fe, contrary to Mg which 



is slightly enhanced compared to Fe. In the LTE abundance analyses, Thevenin (1998) 
did not found any strong enhancement of the alpha elements Ca, Mg, Si, and only a 
slight enhancement for the oxygen element: [O/Fc]~0.1. Then, we shall adopt a moderate 
[(alpha elements + 0)/Fe] value, about 0.1. It results that we adopt for the modeling 
of HD 203608 a surface metallicity of [Z/X]f = —0.50 ± 0.10 with respect to the solar 
metallicity (|Grevesse fe Noels 1993[) . Because this star belongs to the old disk, we shall 
adopt for the calibration an initial abundance Yj=0.25. 

2.3. Mass and radius 

The accurate Hipparcos parallax determination may help to constrain the adopted T e g 
value with the estimate of the radius derived from surface brightness relations, as al- 
ready done for several asteroseismic studies (|Kervella et al. 2 004. Theveni n et al. 2005|) . 
Values of the radius derived from the magnitudes in different infrared bands are given in 



Table [TJ Creevey et al. (2007) have demonstrated the advantage to constrain the stellar 



modeling with asteroseismic data with an independent interferometric measure of the ra- 
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Table 1. Color magnitudes of HD 203608, and inferred radii. 



Band 


Magnitude 


R/ Rq 


V 


4.22 


0.96±0.08 


R 


3.75 


1.06±0.06 


I 


3.45 


1.08±0.05 


J 


3.27 


1.04±0.04 


K 


2.90 


1.08±0.02 


L 


2.85 


1.09±0.02 



Table 2. Physical parameters of HD 203608 



T cB (K) 

[Z/X] f 

L/L Q 

II (mas) 

log g (cms -2 ) 

vsinj (kms -1 ) 



6070±100 
-0.50 ±0.10 

1.39±0.13 
107.98±0.19 

4.30±0.15 
2.4±0.5 
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Fig. 1. Time series of the signal intensity. Night 4 suffered from an unusually large seeing 
(up to 3.5"). The efficiency of the collected signal is highly sensitive to the sky opacity and 
to the seeing. 
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Fig. 2. Time series of the radial velocity of HD 203608 measured with the pipeline reduction 
of HARPS (unfiltered data). The signal is free of any important low- frequency component. 
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Fig. 3. Lomb-Scargle periodogram, and inset of the window function. The time series used 
for this spectrum excludes the noisiest value; night 4 is completely filtered out with such a 
treatment. 



dius. However, the infrared excess and the presence of the disk around the star make this 
estimate inaccurate, as shown for the Vega-like star r Ceti. Its diameter measured with 
interferometry has been first used to constrain its age and mass (|Di Folco et al. 2004[) ; but 



recently, Di Folco et al. (2007) have measured the K band excess due to hot dust exozodi- 
acal disk, and then found a smaller diameter. Therefore, we prefer not to use in this study 
the HD 203608 diameter as a fixed constraint. 



The mass of HD 203608 has been estimated by Jonsell et al. (2005) using evolutionary 
tracks inversion and assuming an age of 14Gyr : M/M© = 0.88. No uncertainty on the 
mass has been proposed by the authors. As for the radius, we cannot consider the stellar 
mass as a constraint. On the other hand, one of the objective of these observations is to 
show the capability of asteroseismology to derive accurate values of these parameters. 



3. Observations 

About 47.6 hours observations were obtained between August 5 and 10, 2006, representing 
2504 individual measurements (Table [3]). The exposure time was 33 s, giving one mea- 
surement each 68.5 s, and a Nyquist frequency about 7.3 mHz, much above the predicted 
oscillation cutoff frequency. For 4 nights out of 5, the observation duration at this winter 
period was greater than 10.3 hours per night. The fourth night was affected by strong 
perturbations, with an extremely bad seeing (up to 3.5") and many cirrus. Both effects 
yielded a low intensity signal (Fig. [l}: most of the photons do not enter the 1" fibre when 
the seeing is too much degraded. 

The mean SNR in the echelle spectrum at 550 nm was typically better than 250 at 
airmass less than 2. The noise level finally derived from the high frequency part of the 
Fourier spectrum is about 5.1 cms" 1 , corresponding to a high frequency noise in the time 
series about 1.6ms" 1 and to a noise equivalent velocity of about 1.0ms" 1 /v / mHz. 

The time series show that HD 203608 does not in fact present any signature of ac- 
tivity (Fig. [2]), contrary to HD 49933 that showed a strong low-frequency modulation 
(|Mosser et al. 2005)) . This difference can be due to the higher stellar age. It can be also 
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Table 3. Journal of radial velocity measurements. The 2nd night was split in two parts. 
AT represents the length of observation each night. The dispersion ctrv is derived from 
the high frequency noise recorded in the power spectrum of each night. 



Date 


Number of 


AT 


ORV 


Aug. 06 


spectra 


(hr) 


(ms- 1 ) 


05 


583 


11.0 


1.5 


06 


358 +65 


6.9 + 1.2 


1.7 + 1.7 


07 


570 


10.9 


2.6 


08 


377 


7.2 


3.0 


09 


551 


10.4 


2.2 



Table 4. Different treatments of the time series, with different filter levels (arbitrary unit, 
similar to the y-axis in Fig. [I]). The dispersion a v is directly measured in the high frequency 
part of the power spectrum; at = <J u \Jn c /2 is the high frequency noise in the time series, 
with n c the number of points in the time series. 



cut 


'/ 


n c 


&t 


a v 


level 


% 




(ms" 1 ) 


(cm s _1 ) 





45 


2054 


2.0 


5.7 


0.1 


42 


2373 


1.8 


5.3 


0.2 


35 


1960 


1.6 


5.1 


0.3 


28 


1580 


1.4 


5.1 


0.4 


21 


1165 


1.4 


5.7 


0.5 


14 


811 


1.3 


6.2 



related to a low value of the stellar inclination, accounting for the low wsinz value. The 
quasi pole-on observations should then hamper any significant signature of the rotational 
modulation. 



4. The seismic signature 

4.1. Power spectra 

Computing the power spectrum of the complete time scries requires a careful treatment 
of the noisiest data. Different ways have been proposed, as the weights introduced in the 



Lomb Scargle analysis by Butler et al. (2004) In fact, this method makes the noise uniform 
in the time series and minimizes the noise level in the power spectrum. As a consequence, it 
modifies the signal and yields a degraded window function, which increases the interferences 
in the power spectrum between the signal, the aliases and the noise. 

On the other hand, eliminating noisy data requires a criterion for defining a threshold, 
and translates immediately in a degradation of the duty cycle. Optimizing simultaneously 
the signal-to-noise ratio (SNR) and the duty cycle is therefore necessary. Choosing the 
optimum threshold for the lowest noise level is possible, but it does not insure the highest 
efficiency for detecting eigenmodes, since highly unpredictable interferences occur between 
the signal, its aliases and the noise. 
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Fig. 4. Echelle diagram, with all peaks with a false alarm probability less than 10%, col- 
lected for 10 different treatments (with increasing cutoff in SNR, hence decreasing value 
of the duty cycle). The size of the symbol is representative of the mode amplitude. A dark 
symbol indicates a peak present in many treatments, whereas a light grey symbol corre- 
sponds to a peak only present in one treatment. The strong peaks at 2.3 mHz with an 
abscissa around 40 /iHz emerge only in the noisiest data set, and therefore are not selected. 
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We have endorsed this impossibility of an a-priori criterion for the best solution, and 
proposed to circumvent it: instead of performing only one single power spectrum, we have 
calculated many, each one corresponding to a different threshold level in the time series 
(Tabled]). Thresholds were based on the signal intensity (Fig. [TJ) . In each case, the LS pe- 
riodogram was computed and the highest peaks were selected according to the analysis in 



Appourchaux (2004) that gives a test for detecting peaks of short-lived p modes embedded 



in a power spectrum of noise. Simulations have confirmed that the single treatment cor- 
responding to the optimum solution is far from providing all significant peaks. The peaks 
with a false alarm probability lower than 10% are plotted in an echelle diagram based on a 
large splitting of 120.5 fiHz (Fig. 2]). The procedure makes then possible the identification 
of the major ridges in the echelle diagram. The most confident peaks in these ridges are 
then identified (Fig. [5]) ; those selected have the minimum false alarm probabilities obtained 
in each periodogram. 



4.2. Asymptotic parameters 

The echelle diagram exhibits clearly the regular pattern constructed by the asymptotic be- 
havior of low-degree high-frequency pressure modes. Clear ridges appear in the frequency 
range [2.3 - 2.8 mHz]. The multiple calculations of periodograms, as many as threshold 
values were fixed, allowed us to put in evidence more peaks than obtained with the cal- 
culation of a single spectrum. For 17 peaks (including the aliases) finally detected in the 
range [2.3 - 2.8 mHz] with a false alarm probability less than 10%, 12 at most were present 
in a single power spectrum. 

The ridges in the echelle-diagram show a curvature in the range [2.3 - 2.8 mHz]. This 
may be interpreted as a significant second order term to the asymptotic law, or as the 
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Fig. 5. Echelle diagram, selecting the major peaks detected in the two major ridges iden- 
tified in Fig. 0] (black symbols), completed with peaks with a larger false alarm probability 
(up to 25%, grey symbols) but following the asymptotic pattern, and excluding the peaks 
present in the neighbor alias of the most prominent ridges. 




-20 20 40 60 80 

Fig. 6. Collapsogram of the echelle diagram, with all identified large spacings rescaled on 
one single value A^ eq , for modes identified in the range [2.2, 2.9 mHz]. Modes £ = or 1 
can be identified, as well as their alias (with a prime). The small spacing derived from the 
frequency difference Si/qi allows us to identify the signature of £ = 2 modes. 



signature of a modulation in the spectrum due to an important density gradient inside the 
star. In order to account for it, we propose a fit of the frequency pattern varying as: 

2 

(1) 



v n ,t = v na +n e Av - £(£ + 1) D 
with 



2 L 



n e = n-n + - 



(2) 



The frequency Do measures the small spacing, as usually done in most of the previous 
works reporting single-site asteroseismic observations. C represents the variation factor 
of the large spacing with the radial order, and defines the local curvature in the echelle 
diagram. It corresponds to a global linear increase of the large spacing with the radial order 
n such as: 



l±.v(n) = A^(no) + C rig 



(3) 
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In order to account for the curvature and the irregularities in the echelle diagram, we 



built a rectified collapsogram, as in Mosser et al. (2008) This rectified collapsogram of the 



echelle diagram (Fig. [5]) , corrected from the variation with frequency of the large spacings 
between modes of same degree, puts in evidence the signature of £=2 modes. Their location 
compared to radial modes is in agreement with the small spacing derived from the frequency 
difference Svqi — v n $ — (is n ,i + ^n-i,i)/2 derived from £=0 and 1 modes. We may then 
propose an identification of eigenmodes with degrees £ = 0, 1 and 2, from which we can 
derive the local values of the large spacing (Fig. [7]) . 

The large spacing is about 120.3±0.5/iHz at 2.6 mHz, matching the value 119 /iHz 



expected from the scaling based on the square root of the mean density y/QM/R 3 . The 
variation factor C, around OA^® ^ fiHz in the range [2.3 - 2.8 mHz], gives a measure of the 
variation of the large spacing with the radial order n. 

4.3. Small spacings 

Variations of the small spacings with frequency show a more complicated pattern than 
expected from simple asymptotics (Fig. [H]) . The decrease with frequency of the small fre- 
quency differences is so pronounced that, instead of the development introduced with Eq[T] 
we prefer a development closer to the asymptotic form (jTassoul 1980|1 : 

£(l+l) A rij „ 

v n ,e = v nQ +n e Au- A Q + ^C (4) 

n + £/ 2 2 

The characteristic frequency Aq introduced in this asymptotic development, around 
24±18/iHz, corresponds at 2.6 mHz to a classical small spacing value Do about 
1.2±0.9/iHz. The correction introduced with the 1/n term appears efficient, but is not 
sufficient to account for the variation of the small frequency difference. 

The evolved star HD 203608 may present a core mainly composed of helium, responsible 
for a significant central decrease of the sound speed. Such a core has a strong influence on 
the second order asymptotic term, through the integral term dc/r, that the classical Tassoul 
development cannot take into account. The rapid density and sound speed variations may 
be modeled as a discontinuity. The asymptotic form in that case (jProvost et al. 19931 see 
their Eq 4 .1) introduces many parameters for describing the modulation of the small 
spacings. With eigenfrequencies identified over a limited frequency range, it has little sense 
to try to fit all of them, so that we prefer to linearize the development in the form: 

£(£+1) n 2 

v n j = v m +n e — — [Pi + Qg n e } + C (5) 

n + £/ 2 2 

It yields small spacings such as: 

Svoi - ~i p i + Qi nt] ~ \c (6) 

5v Q 2 = -[Fi + <?2 nt] (7) 
n 

with Svo2 — Vnfl — ^n-1,2- The parameters introduced by Eq[5]are summarized in Tab. [5) 
Compared to the Aq factor introduced by Eq[H the relative errors on the parameters P± 
and Pi are sensitively reduced. The identification of individual eigenmodes is finally given 
in Table [6] We note that the location of £ = 2 modes is possibly influenced by aliases of 
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Fig. 7. Variation of the large spacing with frequency, for the degrees £ = (squares), 1 
(circles) and 2 (diamonds). The dispersion with respect to the linear fit (dashed line) is 



compatible with the frequency uncertainty around 1.2 /iHz derived from ILibbrecht 1992} 
that yields twice this uncertainty on each frequency difference. 
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Fig. 8. Variation of the small spacings (n + £/2) Svoi/2 and (n+£/2) ^^02/6 with frequency, 
derived from the £ = 1 (circles) and 2 (diamonds) modes, compared to £ = 0. 



the £ = modes, since the configuration of the time series yield reinforced signatures of 
the alias related to a 4-day periodicity. 

4.4. Amplitudes 

According to its mass and luminosity and following the power law given by 
Samadi et al. (2007)| the maximum amplitudes of HD 203608 was supposed to be about 
30 cms -1 . In order to estimate the maximum oscillation amplitude, we have constructed 
synthetic time series, based on a theoretical low degree p-modes eigenfrequency pattern. 
The modes lifetimes are estimated from the eigenfrequency widths, between 1 and 4 /iHz 
FWHM (Hou dek et al. 1999[) . Due to the short duration of the time series, possible large 
uncertainties in the lifetimes estimate translate in very small uncertainty in the result. The 
maximum amplitudes are assumed to follow a Gaussian distribution in frequency. The syn- 
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Table 5. Estimation of the asymptotic parameters, relying on the modes detected in the 
frequency range [2.3, 2.8 mHz], and 3-cr error bars. 



Asymptotic parameters (at 2.6 mHz) 




with a Tassoul-like development 




Av 


120.3±0.5 ^Hz 


C = dAv/dn 


0.4«; 4 3 ^Hz 


A 


24±18 ^Hz 


D ~ Ao/20 


1.2±0.9 mHz 


Second order terms 


Pi 


26.0 ± 7.5 pHz 


P-2 


22.2 ± 3.5 fiRz 


(h 


-5.2 ± 3.2 fiUz 


Q2 


-0.3 ± 1.8 fiHz 


Amplitudes 


^max 


2.6 mHz 


^max 


22±2cms _1 



Table 6. Identified peaks, with V their confidence level (here complementary to the false 
alarm probability), and inferred peaks, as the function of the assumed radial order n'. 
Frequency uncertainty is around 1.2 /iHz, according to the estimated lifetime, the obser- 
vation duration and the SNR (jLibbrecht 1992[) . The detections and identifications outside 
the range [2.3-2.8 mHz] are not as certain as within it. Confidence levels are given when 
modes are directly identified in the echelle spectrum. For £ — 2 modes, they cannot be 
given, since those modes are just inferred assuming a regular echelle pattern. 

1 = L = 1 1 = 2 

n Vobs V ^obs V l^obs 



mHz % mHz % mHz 



17 






2.305 


95 


2.361 


18 


2.369 


91 


2.425 


87 


2.480 


19 


2.488 


79 


2.545 


94 


2.602 


20 


2.608 


86 


2.667 


92 


2.722 


21 


2.728 


94 


2.791 


93 


2.843 


22 


2.849 


88 


2.911 


81 


2.967 


23 


2.973 


71 


3.032 


84 




24 


3.093 


75 


3.154 


71 





thctic time scries arc then calculated using the model of a stochastically excited, damped 
harmonic oscillator (Anderson et al. 1990), and include also a white noise. 

The maximum amplitude is adjusted in order to obtain comparable energy per frequency 
bin in the synthetic and observed spectra (the reference observed spectrum has a cut level at 
0.2, what corresponds to the minimum dispersion a v , as shown in Tabled]). A Monte-Carlo 
approach finally shows that the best agreement is for a signal with a maximum amplitude 
about 22±2 cms -1 rms, with a Gaussian envelope centered at 2.6 mHz and with a 1.2 
mHz FWHM. The simulation shows that the noise component, precisely determined due 
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to the large oversampling of the time series, is 1.57±0.10ms , in agreement with the high 
frequency noise directly determined in the spectrum. 

The observed maximum amplitude and the predicted scaling (33±5 cms -1 rms) agree 
only marginally. The difference may be related to the low metallicity of HD 203608; 
a low metallicity yields a thinner convective envelope, then possibly lower amplitudes. 
Preliminary 3-D simulations of the outer layers of a star with a metallic abundance 10 
times smaller than solar result in a mode excitation rate about 2 times smaller. Actually, 
low metallicity corresponds to a mean opacity smaller compared to the solar one. Then, 
in the super-adiabatic region, where convection is inefficient because of the optically thin 
atmosphere, the radiative flux is larger than in a medium with a solar metallicity. In that 
case, convection can be less vigorous for evacuating the same amount of energy, leading 
to a weaker driving. Therefore, it seems coherent that p modes in a star with sub-solar 
metallicity [Fe/H]~ —0.6 dex present significantly smaller amplitudes than in a star with 
solar metallicity. 

This preliminary analysis has to be refined, since previous asteroseismic subgiant tar- 
gets with much lower or greater metallicity than the Sun did not show such a discrepancy 



compared to Samadi et al. (2007) empirical law {y Ind, [Fe/H]~ — 1.4. ICarrier et al. 20071 
H Ara, [Fe/H]~0.35, IBazot et al. 2005|) . On the contrary, the dwarf target HD 49933 
([Fe/H]~ — 0.37. iMosser et al. 2005]) already showed weaker amplitudes than expected. 



5. Modeling 

Interior models taking into account the new asterosismic constraints were all computed 
using the evolutionary code CESAM2k (jMorel 1997)) . We used the OPAL equation of state 
and the nuclear data from NACRE (Ang ukTet al. 1999] ). The boundary layers are described 
using a model of atmosphere derived from the Kurucz model adapted to an undermctallic 
star (jKurucz 1997[) . For the chemical composition, we used the solar mixture from Grevesse 



& Noels (1993). The revised abundances from Asplund et al. (2005) suggest lower abun- 
dances of C, N, O, Ne and Ar. Many studies showed however that the standard solar models 
using the new abundances of Asplund are in disagreement with the sound speed profile, 
the radius and the helium abundance of the convection zone (jGuzik et al. 2005)) . These 
new abundances need to be confirmed or infirmed by independent 3D NLTE line transfer 
studies of the oxygen element. Nevertheless, we also computed models with the abundances 
of Asplund, and checked that they do not induce any significant change in the results. The 



convection follows the description of Canuto & Mazzitelli (1991) with a mixing length pa- 
rameter A close to 1. A model of the Sun using this description of the convection led to 
A ~ 0.94. This value was then adopted for this study. The stellar models were computed 



with microscopic diffusion using the formalism developed by Burgers (1969) 



5.1. Description of the fit 

The models are constrained using the following observational quantities : Toff, L/L@, 
[Z/X]f, and asteroseismic global parameters derived from the asymptotic development. 
Temperature, luminosity and final metallicity are fixed in the code; all other parameters 
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are free, including the mass and the age. The difference between the computed models and 
the observations is quantified by the x 2 function: 



X 



s /,„obs „mod^ 2 
2 / Pi P 



i=l 



Best-fit models are found by minimizing this function. We then obtain an estimate of the 
parameters which cannot be measured by observations, such as the mass of the star, its 
age, the initial abundance of helium Y and the initial metallicity [Z/X]i. 

The value of the initial helium content cannot be measured, but it can be assessed using 
the relative helium to metal enrichment of the galaxy AY/AZ through the relation: 
AY 

Y =AZ Z + Y * ® 



where Y p is the primordial helium content. Pagel & Portinari (1998) have shown that the 
enrichment was such that 2 < AY/AZ < 5. With a primordial helium abundance of 
Y p = 0.238, the initial helium content in the star is found between 0.25 and 0.27. 



5.2. Results 

5.2.1. Taking into account the mean large spacing 

The observed frequencies from the data analysis range from 2.3 mHz to 3.2 mHz. The large 
spacing derived from the observations is determined with an uncertainty of 0.5 /iHz in this 
frequency range, for the modes I = 0, 1 and 2. The first models had to fit the values of T e g, 
L, [Z/X]f, [Z/X]i and Yj, plus this single asteroseismic constraint: only models which fit the 
mean large spacing with an accuracy better than 0.5 /xHz were kept. The resulting position 
of the star in the HR diagram is narrowed (Fig. [5]) and the precision of the parameters of 
the star is improved (Table [7]). 



5.2.2. Taking into account the second order terms 

Taking into account the other asteroseismic parameters is required to constrain the models 
more efficiently We therefore use the development expressed by Eq which allows us 
to take into account the variation of the large spacing with frequency, and the different 
behaviors of the small spacings Svoi and 61/02 (Eqs[5]and[7]). The observational asteroseismic 
constraints for the stellar modeling are then: Av : C, P±, Q\, P2 and Q2- 

Figure [TU] shows for instance the evolution of the large spacing with frequency, both 
for the observations and for one of the best-fit models. These evolutions are apparently in 
very good agreement: the values of the large spacing at 2.6 mHz agree within 0.1% and 
the slope of the variations of the large spacing within 10%. The only parameter which 
shows a marginal agreement is the one describing the evolution of Svqi, i.e. Q\, as can 
be seen on Fig. [TTJ In fact, none of the models reaches the value of the slope derived 
from the observations. The agreement in the case 2-0 simply derives from the fact that 
the perturbation due to a dense core depends mainly on the factor n + £/2, and does not 
affect significantly the spacing 8 v$2- On the other hand, it strongly affects the spacing Svqi- 
Hence, a small discrepancy between the best fit model and the observation translates into 
a large discrepancy on this small spacing. 
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Fig. 9. Position in the HR diagram of the models which fit the mean value of the large 
spacing. The inner frame represents the 1-er limits of log(T c ff) and log(L/LQ). The dashed 
lines set the 1-cr limits of the large splitting. The dotted line represents the location in the 
HR diagram of the models having a mean value of <5^o2 = 6.5/iHz, like the one derived from 
the observations. The solid line shows an example of an evolutionary track fitting Ai>. Its 
characteristics are: M = 0.94 M , age = 6.9 Gyr, Y x = 0.26, [Z/X] s = -0.38. The numbers 
along the evolutionary track stand for the age of the star in Gyr. 
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Fig. 10. Variation of the large spacing with frequency, for the £=0 (squares), 1 (circles) 
and 2 (diamonds) modes. The full symbols and the dashed line stand for the observations. 
The open symbols and the solid line stand for a model fitting the observational constraints. 



This phenomenon mainly accounts for the impossibility to get very low values for x 2 
(the minimum value is Xmin — 9-5)- These observational constraints yield an estimate of 
the stellar parameters, which are given in the second part of table [71 Fig. rT3] shows the 
position of HD 203608 in the HR diagram resulting from the minimization. 

Even if sensitive diffusion effects are expected for such a target, models were also com- 
puted without adding microscopic diffusivity, in order to study the impact on the evolution 
of the star. Agreement is worse, the best model providing Xmin — 12-3. The resulting mod- 
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Fig. 11. Variation of the second order terms with frequency, derived from the I — 1 (circles 
and dashed lines) and 2 (diamonds and dotted lines) modes, compared to I = 0. The open 
symbols represent the model and the full symbols, the observations. Agreement is better 
in the case 2-0 than in the case 1-0, since the perturbation due to a dense core depends 
mainly on the factor n + 1/2. 

Table 7. Physical and seismic parameters of HD 203608 derived from the modeling of the 
star (all models calculated with the convection parameter A = 0.94). The models of type 
a are those which only take into account the mean value of the large spacing. The models 
of type b are those which use the fitted values of Av, C, [Z/X]i, Pi, Qi, P2 and Q 2 - 



Physical parameters models a models b 

T eff 6050 ± 100 K 6037 ± 19 K 

L/Lq 1.39 ± 0.13 1.341 ± 0.016 

[Z/X] f -0.55 ± 0.05 dex -0.54 ± 0.03 dex 

Yi 0.26 ± 0.01 0.252 ± 0.003 

[Z/X]i -0.42 ± 0.05 dex -0.416 ± 0.03 dex 

M/Mq 0.93 ± 0.07 0.939 ± 0.009 

age 7.5 ± 2.6 Gyr 7.25 ± 0.07 Gyr 

R/Rq 1.064 ± 0.014 b mean value over the frequency 

Seismic parameters models a h models fc" 

Az/ 120.3 ± 0.5 /jMz 120.4 ± 0.4 /jMz 

8v 02 6.9 ±1.5/iHz 6.75 ± 0.05 /jMz 

C 0.290 ± 0.003 fiRz 

Pi 32.54 ± 0.48 fiHz 

P 2 22.50 ±0. 14 fiRz 

Qj -1.78 ± 0.06 fiHz 

Q 2 -0.33 ± 0.05 fiRz 
range [2.3 - 3.2 mHz] 

» at 2.6 mHz 



els are expectedly older than those including diffusion (about 8 Gyr old), and the effective 
temperature is about 100 K higher. 
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Fig. 12. Minimization of the x 2 function. The minimum value of x 2 is Xmin — 9.5. The 
solid line, the dashed line and the dotted line delimit the areas in the HR diagram at 1, 2 
or 3 a where respectively \ 2 < Xmin + 1,4 and 9. As in Fig. [9j the inner frame represents 
the 1-ct uncertainties previous to this asteroseismic run. 

Finally, our study yields to a star reaching the end of the main sequence. The hydrogen 
is almost entirely exhausted in the center: the mass fraction of remaining hydrogen in 
the core is of about 13%. The best models show no trace of a convective core during 
the evolution sequence, except plausibly during the first billion of years. As explicited by 
Tables [5] and the precision on the stellar parameters has been significantly improved. 
Localization in the HR diagram is refined by a factor greater than 5 in temperature and 
about 8 in luminosity. The error bar on the mass is now defined, as low as 0.01 M & . The 
precision on the age, is much better, and we note that the age of the star corresponds to 
the low limit given by previous works (6.5— >14.5 Gyr according to |Bryden et al. 2006| ). 

6. Conclusion 

This single-site 5-day long asteroseismic run on HD 203608 has given a much more precise 
view of this star of the old galactic disk. We have developed a method for extracting the 
peaks with the lowest false alarm probability. This method proves to be efficient for single- 
site observation with rapidly varying photometric conditions. The performing of multiple 
periodograms combined with a statistical test allows us to extract more peaks than with a 
single treatment. Using a criterion combining the minimization of false alarm probabilities 
and echelle diagram analysis, we have identified 15 t = and £ — 1 eigenmodes in the 
spectrum of HD 203608, from which we have derived lower amplitudes modes (including 
£ — 2 modes) as well as the asymptotic parameters. 

Despite the very short duration of the run, yielding a limited precision for the identi- 
fied eigenfrequencies, the fitting of the spectrum has required a more precise development 
than the usual second order term — £(£ + 1) Do. We have shown that the Tassoul original 
form with a second order term decreasing with frequency — £{l + l)/(ra + £/2) Aq must 
be preferred to the development —£(£ + 1) D often used for interpreting ground-based 
observations. Nevertheless, the asymptotic development cannot account for the precise os- 
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dilation pattern of HD 203608: this star exhibits clearly a large spacing with a sensitive 
dependence on the radial order, and small spacings depending on the mode degree. This 
dependence observed in the data are confirmed in the modeling: the strong composition 
and sound speed gradients in the small core mainly composed of helium are responsible of 
the modulation of the oscillation pattern. The eigenfrequency precision in our data set is 
not accurate enough to give additional independent constraints on the core size; continuous 
long-duration observations are required for such a task. 

Contrary to a similar F type dwarf target (HD 49933. [Mosser et al. 2005]) . HD 203608 
does not exhibit any noticeable activity. This may be due to the geometric configuration 
of the observation, with a possible very low inclination axis. By now, HD 203608 presents 
the lowest metallicity among dwarfs observed in asteroseismology. Similarly to HD 49933, 
modes amplitudes are sensitively smaller than expected from the scaling law. Two effects 
may explain this: first, both stars are undermetallic; second, the scaling has not yet been 
calibrated on dwarf F stars. Observations with the satellite CoRoT will help understanding 
that behavior. 

The modeling of HD 203608 has been achieved with the evolution code CESAM2k 
(jMorel 1997P - Taking into account the asteroseismic constraints (large and small spacings) 
allows us to propose a much more precise description of this star. Error bars on the phys- 
ical parameters have been divided by a factor of 2 to 8, in the framework of the present 
physical description used in this work. The age we derive for HD 203608 is about 7.25±0.07 
Gyr. Hydrogen is almost exhausted in the core. With improved values of L, T e g and log g 
(4.356±0.016), a better signature of the metallic abundances is certainly possible and may 
help to improve the localization of this star in the thin or thick galactic disk. 
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